Global Attractivity of the Equilibrium 
of a Difference Equation: 

An Elementary Proof Assisted by Computer Algebra System 



oo 
O 
O 
(N 

O 

Q 



Orlando Merino 
Department of Mathematics 
University of Rhode Island 
Kingston, RI 02881 
merino @math. uri. edu 



Abstract 



Let p and q be arbitrary positive numbers. It is shown that if q < p, then all solutions to the 
difference equation 



x n+l = P + qXn n = 0,1,2,..., x_i>0, x >0 (E) 
g 1 + X n -l 

converge to the positive equilibrium x = \ I q — 1 + y/(q — l) 2 + 4p J . 



parameters, thus completing the proof of a conjecture of Ladas. 



The above result, taken together with the 1993 result of Kocic and Ladas for equation (E) 
ON with q> p, gives global attractivity of the positive equilibrium of (E) for all positive values of the 
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1 Introduction and Statement of Main Theorem 

The difference equation 

x n+ i = f qXn , n = 0,1,2,..., x-i > 0, x Q > (E) 

has a unique equilibrium x = | — 1 + \/{q — l) 2 + 4pJ . In 1993, V. Kocic and G. Ladas proved 

in [3] (namely, part (e) of Theorem 3.4.3) that for q > p > 0, all solutions to Eqn.(E) converge to 
the equilibrium. Since then the region of positive parameters p, q for which global attractivity of 
the equilibrium holds has been extended in several publications: V. Kocic, G. Ladas, and I. W. 
Rodrigues [I], C. H. Ou, H. S. Tang and W. Luo [T3J, H. A. El-Morshedy [TJ, R. Nussbaum [H], 
V. Jimenez Lopez [2]. However, there is a large region of parameters that is not covered by the 
results obtained in those publications (see Figure [T]). Indeed, a well known conjecture of Ladas 
states the positive equilibrium to (E) is a global attractor for all positive values of the parameters 
p and q, (for example, see [1UJ . [BJ, [7]). The authors in reference [TT] claim to have proved the 
conjecture, but there is a mistake in their proof, see [2] for details on this. 

The main result of this work, given below, establishes global attractivity of the equilibrium of 
(E) for q < p, thus completing the proof of the conjecture. 

Main Theorem // < q < p, then all solutions to the difference equation (E) converge to the 
positive equilibrium x. 




Figure 1: Boundaries of regions in parameter space. Key: (a) q > p (Kocic-Ladas, 1993), (b) 2(g+l) > p 
(Kocic-Ladas-Rodrigues, 1993), (c) 2(q 3 - q 2 + q + ^g 4 - 1 - l)/(q - l) 2 > p (Ou-Tang-Luo, 2000), 

(d) \ (q + y/(q - l) 2 + Ap - l \ < (El-Morshedy, 2003, and Nussbaum, 2007) (e) Ap (q - l) 2 > 25 

Jimenez Lopez (2008). It is interesting to note that Nussbaum obtained El-Morshedy' s bound expressed 
in different form. Also, Jimenez Lopez' region is an improvement on regions obtained previously by 
Ou-Tang-Luo, El-Morshedy, and Nussbaum for q near 1 and p > 112.1 (not shown in the figure). The 
point (p, q) = (20, 4) is in the exterior of all regions bounded by the different curves. 



The study of Eq.(E) is facilitated by the introduction of a suitable change of variables, and 
thus we begin with it. The substitutions 

a 1 
V= T^ q= A 

in Eq.(E) yield the equation 

a + y n 



and xi = qye (£ = —1,0,2,...) 



(1) 



Vn+l 



A + y n -i 

The equilibria x of Eq.(E) and y of Eq.Q satisfy 

x = qy 

One may view Lyness' equation (see [3]), 

a + z n 



n = 0,1,2,..., y—i, yo € (0, oo) 



Zn+l 



Zn-1 



, n = 0, 1,2, . . . , z e (0,oo) 



(2) 



(3) 



(4) 



as a limit case of ^ that results by setting A = 0. Lyness' equation is known to posess an invariant 
function i.e., a continuous, nonconstant real valued function g(x,y) defined on (0,oo) x (0,oo) 
such that every solution {z n } to Q satisfies 

g(z n+ i,z n ) = g(z ,z-i), n = 0,1,2,... (5) 

A formula for g(x,y) was found by Lyness [12] for 5 = 1, and for positive a by Ladas [9]: 

(l + x)(l + y)(a + x + y) 



g(x,y) 



xy 



(6) 



It is easy matter to show that g(x, y) has a unique critical point on the positive quadrant, namely 
(z, z), where z is the equilibrium of Eq.Q. Since g(x, y) is large near the boundary of the positive 
quadrant, g(x,y) has a strict global minimum which is is attained at (z, z) (see [16], [5]), i.e., 



g(z, z) < g(x, y) for (x, y) G (0, oo) x (0, oo), (x, y) / (z, z). 
The function g(x,y) plays a fundamental role in our proof. See Figure [2} 



(7) 





Figure 2: 3D-Plot and contour plot of the invariant function g(x,y) for a = 2. 

In Section[2]we prove several lemmas before giving the proof of the main theorem, which is done 
in order to simplify the exposition. For practical reasons, some of the calculations were performed 
with the computer algebra system (or CAS) Mathematica |15| . Code for such calculations, written 
in the Mathematica language, is given in Section [4] 
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2 Proofs 

For typographical convenience we will use the symbol u to represent the equilibrium y of Eq. ^ . 
By direct substitution of the equilibrium u = y into Q we obtain 

a = u 2 + {A-l)u (8) 

By Q, a > implies u > 1 — A. Using ([8]) to eliminate a from ([2| gives the following equation, 
equivalent to (E): 

u 2 + (A — 1) u + y n 

2M+1 = 1— , ra = 0,1,2,..., y_i, yo G (0, oo) (9) 

-A + 2/n-l 

Therefore, to prove the main theorem it suffices to prove that all solutions of Eq.([9| converge to 
the equilibrium u, and this is what the rest of the proof is geared to do. 
The following statement is crucial for the proof of the main theorem. 

Lemma 1 u > 1 if and only if q < p. 

Proof. Since x = qy = qu, we have u > 1 if and only if x > q, which holds if and only if 
^ (q — 1 + sj {q — l) 2 + 4pj > q. After an elementary simplification, the latter inequality can be 
rewritten as q < p. □ 

The hypothesis of the main theorem states q < p. In view of Lemma [T] we assume in the rest 
of the exposition that 

u > 1 (10) 

For the equilibrium u of Eq.Q to be also the the positive equilibrium of Lyness' equation Q one 
must have 

a + u 

u = — (11) 

that is, 

a = u 2 -u (12) 
For the value of a found in (12) the invariant function ^ becomes 

= (l + x){l + y){u 2 -u + X + y) 
xy 

One would hope that g(x, y) turns out to be a Lyapunov function for Eq.Q (see j5]). We shall see 
in the proof of Lemma [2] that this is not the case. Nevertheless we will be able to use the function 
g(x,y) to complete the proof of the main theorem. First we need some elementary properties of 
the sublevel sets 

5(c) :={(*,*) :g(s,t) <c}, c>0 
We denote with Qe(u, u), £ = 1, 2, 3, 4 the four regions 

Ql(u,u) := {(x,y) G (0,oo) x (0, oo) : u < x, u < y } 

Q2(u, u) := {(x, y) G (0, oo) x (0, oo) : x < u, u < y } 

Q 3 (u,u) := {{x,y) G (0, oo) x (0, oo) : x < u, y < u } 

Qa(u, u) := {(x, y) G (0, oo) x (0, oo) : u < x, y < u } 
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Let 



T(x,y) 



y, 



u 2 + {A-l)u + y 



(x,y) G (0,oo) x (0,oo) 



A + x 

be the map associated to Eq. ([9]) (see [S]). 

Lemma 2 // (x, y) G Q2(u, u) U Qi(u, u) \ {(u, u)}, then g(T(x, y)) < g(x, y). 
Proof. Set 



A calculation yields 
Ai(x,y) 



Ai(x, y) := y) - g(T(x, y)) 

A (u — u 2 + ux — y) (u — Au — u 2 + Ax + x 2 — y) (1 + y) 
x (A + x) y (—u + Au + u 2 + y) 



(14) 



(15) 



(16) 



By inspecting the factors in (16 ) one can see that Ai(x, y) changes sign on the line y = u — u+ux 
and on the parabola y = u — Au — u 2 + Ax + x 2 . Both curves are the graphs of strictly increasing 
functions of x that intersect in the positive quadrant {(x,y) : x > 0, y > 0} at a unique point, 
namely (u, u). At the point (u,u), the slope of the parabola, 2u + A, is larger than the slope of 
the line, u. See Figure[3} Clearly Ai(x,y) becomes negative for points (x,y) between both curves, 
which, except for (n, u), are contained in the complement of Q2(u, u) U Q±(u, u). □ 




Figure 3: Ai(x, y) < for (x, y) in the interior of the shaded region bounded by the line y = u — u 2 + ux 
and the parabola y = u — Au — u 2 + Ax + x 2 . 



Lemma 3 // (x, y) G Qi(u, u) U Q 3 (u, u) \ {(u, u)}, then g(T 2 (x, y)) < g(x, y). 
Proof. Set 



A 2 (x,y) := g(x,y) - g{T 2 (x,y)) 



(17) 



A calculation yields 



-A 2 u 3 + A 3 u 3 + A 4 u 3 - A 5 u 3 H 386 monomials \- 2 A 2 u x 2 y 4 " + Au x 3 y 4 + Ay 5 

x (A + x) y (A + y) (—u + A u + u 2 + y) (—u + A 2 u + u 2 + Au 2 — u x + Au x + u 2 x + y) 

(18) 



Since u — u > by (10), the denominator of A 2 (x, y) is positive. We now proceed to prove that 



the numerator of ^{x, y), denoted by Num, is positive on Qi(u, u) \ {(u, u)}. Begin by changing 
variables: x = xq + u, y = yo + u, where xo > 0, yo > to get 

Num = A n 4 x 2 + A u 5 x\ + 2 A u 6 Xq + 2 A u 7 Xq + 2 A u 3 Xq - A u A x|] H 287 monomials (19) 



Since some monomials in the right-hand-side of (|19j) have negative coefficient, it is not obvious 
that Num > for xq > 0, yo > 0, (xo,yo) 7^ (0,0). The last step here consists in considering 



subcases (yo > xq, yo < xq, and yo = xq). Substituting yo = xq + k and u = 1 + t for k > and 
t > 0, we obtain 

Num = 2 Ak 2 + 8 A 2 k 2 -\ 368 monomials \-2Atxl (20) 

By using a computer algebra system one can show that the coefficients of the monomials that form 



the expression in the right-hand-side of (20) are all positive integers. Thus Num > whenever (a) 
xq > 0, yo > xq. It can be easily verified that Num > in the other cases: (b) yo > 0, yo < xq, 
(c) xo = 0, yo > 0, and (d) xo > 0, yo = 0. We conclude that A2(x,y) > on Qi(u, u) \ {(u,u)}. 
See Section U] for details. 

To verify that Num > on Qs(u, u)\{(u, it)}, we consider first Num on the interior of Qa(u, u), 
and introduce the algebraic transformation 

uw uv . . . . . . 

y = — -, (w,v) € (0,oo) x (0,oo) 



W + 1 V + 1 

that maps the open positive quadrant {(w, v) : w > 0, v > 0} one-to-one and onto the interior of 
Qs(u,u). Then substitute n = 1 + t with t > 0, and consider subcases above, below, and on the 
diagonal just as we did before to conclude by direct inspection of coefficients that the transformed 
expression for Num is positive. 

A similar strategy works for proving A2(x,y) > on the open line segment with endpoints 
(u, 0) and (it, it), and on the open line segment with endpoints (0,it) and (it, it). See Section [4] for 
details. □ 

Proof of the Main Theorem. Let (</>, tp) £ (0, do) x (0, do). Let {y n }n>-i be the solution to 
([9]) with initial condition (y_i,yo) = (^VO) an d let {T n (</>, ip)} n >o be the corresponding orbit of 
T. Define 

c:=liminf 5 (T n (0,VO) (21) 

n 

Note that c < 00, which can be shown by applying Lemma [2] and Lemma [3] repeatedly as needed 
to obtain a nonincreasing subsequence of {g(T n (4>, tp))} n >o that is bounded below by 0. Let 
{g(T nk (4>, V0)}fc>o be a subsequence convergent to c. Therefore there exists c > such that 

g(T nk ((j),ip)) < c for all A; > 0, 

that is, 

T nk ((f), if;)) e 5(c) := {(s,t) : 5 (s,t) < c} for k > 
The set 5(c) is closed by continuity of g(x, y). Boundedness of 5(c) follows from 

n . (1 +x)(l +y)(u 2 -u + x + y) , . . 

< x, y < =g(x,y) = c, for (x,y)e5(c) 

xy 

Thus 5(c) is compact, and there exists a convergent subsequence {T nk e(cj),ip))}£ with limit (x,y) 
(say). Note that 

C = lim y(rM<^)) =gM) (22) 

£^00 

We claim that (x,y) = (u, it). If not, then by Lemma [2] and Lemma |3j 

min{y(T(x, y)),g{T 2 {x, y))} < c (23) 



Let || • || denotes the euclidean norm. By (23) and continuity, there exists 5 > such that 

\\(s,t) - (x,y)\\ < 5 =^ mm{g(T(s,t)),g(T 2 ( S ,t))} < c (24) 
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Choose L € N large enough so that 

||r^(<£,v>)-(£,y)|| <6 (25) 



But then (|24j) and (|25j) imply 

min{< 7 (T"^+ 1 (0,^)),5(T n ^ +2 ( S ,t))} < c (26) 



which contradicts the definition (21 ) of c. We conclude (x, y) = (u, u). From this and the definition 
of convergence of sequences we have that for every e > there exists L E N such that ||T n *£ ((ft, tft) — 
(u,u)\\ < e. Finally, since 

maxj|y nfcL _i - u\, \y nkL -«|} < ||(Vn* £ -l ~ «, Vn kh - u)\\ = \\T nh L - ( u , u)\\ 

we have that for every e > there exists L £ N such that \y nkL — u\ < e and \y Hk _i — u| < e. 
Since u is a locally asymptotically stable equilibrium for Eq.Q (this is Lemma 3.4.1 in page 67 
of [3]), it follows that y n — ► u. This completes the proof of the Main Theorem. □ 

Acknowledgment The author thanks M. Kulenovic for helpful discussions and for suggesting 
the change of variables ([T]) which established a direct connection of the main difference equation 
to Lyness' equation, and J. Montgomery for suggesting a simplification of the proof. 
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4 Appendix: Mathematica Code 
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ln[l]:= g[{x_, y_}] = 



(1 + x) (1 + y) (a + x + y) 



/ . {a -» u 2 - u} ; 



xy 



ln[2]:= T[{x_, y_} ] = {y, 



u 2 + (A - 1) u + y 



}» 



A + X 



ln[3]:= DELTA2 = g[{x, y} ] -g[T[T[{x, y}]]]; 

ln[4]:= Num = Numerator [ Together [ DELTA2 ] ] ; 

ln[5] := Numll = Expand [Num / . {x-»xO + u, y -> yO + u} ] ; 

ln[ 6] := Min [Expand [Numll / . {yO-»xO + k, u -» 1 + t} ] /. Plus -» List / . {xO -» 1, k -> 1, t-»l,A-»l}] 

ln[ 7] := Min [Expand [Numll / . {xO-»yO + k, u -» 1 + t}] /. Plus -» List / . {yO -» 1, k -> 1, t-»l,A-»l}] 

ln[ 8] := Min [Expand [Numll / . {xO -» yO, u -» 1 + t} ] / . Plus -» List / . {yO -» 1, t -» 1, A -» 1} ] 

ln[9] := Min [Expand [Numll / . {xO -» 0, u -» 1 + t} ] /. Plus -» List / . {yO -» 1 , t -» 1 , A -> 1} ] 

in [10]:= Min [Expand [Numll / . {yO -> 0, u -» 1 + t}] /. Plus -> List / . {xO -» 1, t -» 1, A -> 1} ] 

ln[ll] := Num22 = Numerator [Together [DELTA2 / . {x -> u w / (w+1), y-tuv / (v + 1) }] ] ; 

in [12] := Min [Expand [Num22 / . {u -> 1 + t, v -> w + k} ] /. Plus -» List / . {k -> 1, w -» 1 , t -» 1, A -> 1} ] 

ln[13]:= Min [Expand [Num22 / . {u-»l+t, w-»v + k}] /. Plus -» List / . {k -> 1, v -> 1 , t -> 1 , A -» 1} ] 

ln[ 14 ] := Min [Expand [Num22 / . {u -> 1 + t, v -» w} ] / . Plus -» List / . {w -» 1, t -» 1, A -» 1} ] 

in [15] := Num22a = Numerator [Together [DELTA2 / . {x -> u w / (w+1), y -» u} ] ] ; 

ln[16]:= Min [Nura22a / . Plus -> List / . {u -> 1 , A -» 1, w -> 1} ] 

ln[17]:= Num22b = Numerator [Together [DELTA2 / . (y-tuv / (v+1), x-»u}]]; 

in [18]:= Min [Expand [Num22b / . u -> 1 + t] / . Plus -> List / . {u -> 1, A -> 1 , v -> 1 , t -» 1} ] 

Table 1: Mathematica code needed to do the calculations in Lemma |3j This input was tested on 
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Input # 



Description of input 



In[l] 
In [2] 
In [3] 
In [4] 
In [5] 
In [6] 
In [7] 
In [8] 
In [9] 
In [10] 
In[ll] 
In [12] 
In [13] 
In [14] 
In [15] 
In [16] 
In [17] 
In [18] 



Define the invariant function for Lyness' equation. 
Define the map associated to difference equation (top. 



Define the expression A 2 (x,y) given in (17). 



Define word Num to be numerator of expression A 2 (x, y). 

Replace x and y in Num by x$ + u and y^ + u. 

Calculate smallest coefficient of monomials in Numll for yo > xq. 

Calculate smallest coefficient of monomials in Numll for xq > yo- 

Calculate smallest coefficient of monomials in Numll for y = xq. 

Calculate smallest coefficient of monomials in Numll on line {(0,?/o) : Ho > 0}. 

Calculate smallest coefficient of monomials in Numll on line {(xq,0) : Xq > 0}. 

Numerator of map that takes mt Qz(u,u) to the positive quadrant in the plane. 

Calculate smallest coefficient of monomials in Numll for v > u. 

Calculate smallest coefficient of monomials in Numll for v < u. 

Calculate smallest coefficient of monomials in Numll for v = u. 

Numerator of map that takes (open) line segment joining (0, u), (u,u) to (0, oo) 

Calculate smallest coefficient of monomials in Numlla. 

Numerator of map that takes (open) line segment joining (u, 0), (u,u) to (0, oo) 
Calculate smallest coefficient of monomials in Num22b. 



Table 2: Explanation for Mathematica input that appears in Table [TJ 



